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The Optical Gravitational Lensing Experiment. Catalog of 
stellar proper motions in the OGLE-II Galactic bulge fields 
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ABSTRACT 

We present a proper motion (/i) catalogue of 5,080,236 stars in 49 Optical Gravitational 
Lensing Experiment II (OGLE-II) Galactic bulge (GB) fields, covering a range of 
— 11° < I < 11° and —6° < b < 3°, the total area close to 11 square degrees. The 
proper motion measurements are based on 138 — 555 /-band images taken during 
four observing seasons: 1997-2000. The catalogue stars are in the magnitude range 
11 < I < 18 mag. In particular, the catalogue includes Red Clump Giants (RCGs) 
and Red Giants in the GB, and main sequence stars in the Galactic disc. The proper 
motions up to /i = 500 masyr _1 were measured with the mean accuracy of 0.8 ~ 3.5 
masyr -1 , depending on the brightness of a star. This catalogue may be useful for 
studying the kinematic of stars in the GB and the Galactic disk. 



Key words: Galaxy:bulge - Galaxyicenter 
Galaxy:structure - astrometry 



Galaxy:kinematics and dynamics- 



1 INTRODUCTION 

The Galactic Bulge (GB) is the nearest bulge in which in- 
dividual stars can be studied in detail. A study of stellar 
populations and stellar dynamics in the Bulge may help us 
understand how bulges formed, what are their populations, 
gravitational potential and structure. 

The proper motions with precise photometry might 
make it possible to separate the observed populations based 
on their kinematics. Such study has first been done by 
ISpaenhauer. Jones fc Whitford I jl992) with photographic 
plates only for a few hundred of brightest red giants in 
Ba ade's window. R e cently deeper study has been done 
by iKuiiken fc Rich! i2002l) with Hubble Space Telescope 
(HST)/WFPC2 in Baade's window. 

Several groups have carried out gravitational microlens- 
ing observations toward dense stellar fields, such as the Mag- 
ellanic clouds, the Galactic center and disc. Until now, hun- 



dreds of event s have been f o und ( EROS: lAubour g et all 
119931: OG LE: lUdalski et alJ l200d : IWozniak et alJ 12001 1 
MACHO: [Alcock et alJ l2000t MOA: iBond et alJ 1200 it 
ISumi et al.ll2005t . and thousands are expected in the up- 



coming years by MOA 1 , OGLE-III 2 and other collabora- 
tions. 

It is well known that the gravitational microlensing 
survey data is well suited for numerous other scientific 
projects fsee |Paczvnsk]ll99d : lGouldlll996l) . The studies of 
the Galactic structure certainly benefit from this type of 
data. The microlensing optical depth probes the mass den- 
sity of compact objects along the line of sight and the 
event time-scale distribution is related to the mass func- 
tion and kinematics of the lensing objects. Observed high 
optical depth may be explained by the presence of the bar 
JUdalski et alJ 11994 lAlcock et all [l997 . 120001 ISumi et alJ 
120031: 1 Afonso et al.ll2003t iPopowski et alll2003l) There is sub- 
stantial evidence that the Galaxy has a bar at its center 
ide Vaucouleurs! 11964 iBlitz fc Spergel Il99lt IStanek et alJ 
11994 Il997t iKiraea fc Paczvriski 11994 lHafner et al.ll200(tl . 
However, the parameters of the bar, e.g., its mass, size, and 
the motion of star s within it, still remain poorly constrained. 

IStanek et alJ jl997f) used the Red Clump Giants 
(RCGs) to constrain the axial ratios and orientation of 



1 http : //www.phys . canterbury . ac .nz/~physib/alert/alert . html 



http : //www. astrouw . edu. pl/"ogle/ogle3/ews/ews .html 
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the Galactic bar. These stars are the equivalent of the 
horizontal branch stars for a metal- rich population, i.e., 
relatively low-mass core helium burning stars. RCGs in 
the GB occup y a distinct regio n in the colour magni- 
tude diagram ( Sta nek et alJ 12000 and references therein). 
The intrinsic width of the luminosity d istribution of RCGs 
in the GB is small, abo ut 0.2 mag dStanek et al.l Il997t 
iPaczvriski fc StanekHl99^ . Their observed peak and width 
of the luminosity function are related to the distance and 
radial depth of th e bar. 

Furthermore, iMao fc PaczvnskH (I2002T) suggested that 
there should be a difference in average proper motions of 1.6 
masyr -1 between the bright and faint RCG sub-samples, 
which are on average on the near and the far side of the 
bar, respectively, if their tangential streaming motion is 100 
kms - 1 . Following this suggestion, ISumi. Ever fc Wozniakl 
( 2003) measured mean proper motion of bright and faint 
RCGs in one OGLE-II field in Baade's Window, and they 
found a difference to be 1.5 ± 0.11 masyr -1 . 

To expand this analysis we measured proper motions in 
all 49 GB fields observed by th e Optical Gravitatio nal Lens- 
ing Experiment 3 II ( OGLE-II; Hjdalski et all2000ft for stars 
down to / = 18 mag, which is sufficiently deep to include 
RCGs. There are several earlie r proper motion cat alogues of 
this ge neral area (c.f. USNO- B lMonet et alJl200sl Improve d 
NLTT lSalim fc Gould lEool and Tvcho-2 lHe)g et aljEoobl . 
Though the area covered by our catalogue is relatively small, 
it reaches deeper and covers a wide range of proper mo- 
tion (fi < 500 masyr -1 ) with the accuracy as good as (~ 
1 masyr -1 ). In §|3| we describe the data. We present the 
analysis method in §|3]and0] In §|SJ|n]and|7]we describe the 
property, zero point and problems in our catalogue. Discus- 
sion and conclusion are given in §0 



2 DATA 

We use the data collected during the second phase of the 
OGLE experiment, between 1997 and 2000. All observations 
were made with the 1.3-m Warsaw telescope located at the 
Las Campanas Observatory, Chile, which is operated by the 
Carnegie Institution of Washington. The "first generation" 
camera has a SITe 2048 x 2048 pixel CCD detector with 
pixel size of 24 fim resulting in 0.417 arcsec/pixel scale. 
Images of the GB fields were taken in drift-scan mode at 
"medium" readout speed with the gain 7.1 e - /ADU and 
readout noise of 6.3 e - . A single 2048 x 8192 pixel frame 
covers an area of 0.24 x 0.95 deg 2 . Saturation level is about 
55,000 ADU. Details of the instrumentation setup can be 
seen m lUdalski. Kubiak fc Szvmariski I Jl997f) . 

In this paper we use 138-555 /-band frames of the 
BUL_SCl-49 fields. The centers of these fields are listed in 
TableQ The time baseline is almost 4 years. There are gaps 
between the observing seasons when the GB cannot be ob- 
served from the Earth, each about 3 months long. The me- 
dian seeing is ~ 1.3". We use t he VI photometric maps of 
the standard OGLE template dUdalski et al.l I2002T) as the 
astrometric and photometric references. 

3 see [http://www.astrouw.e du.pl/~oglel or 
http : //bulge .princeton. edu/ "ogle 



Only about 70% of the area of the BUL.SC1 field over- 
laps with the extinction map made bv IStanekl |l996). The 
extinction map covering all OGLE-II fields has been con- 
structed bv lSumil J2003T . 



3 ANALYSIS 

The a nalysis in this work follows ISumi. Ever fc Wozniakl 
( 2003), except our procedure makes it possible to detect high 
proper motions, extends to the limiting magnitude down to 
/ = 18 mag and corrects the s ystematic effects. Th e stan- 
dard OGLE template given bv lUdalski et, al] l)2002f) serves 
as the fixed astrometric reference in our analysis. In order to 
treat properly frame distortions in the y-axis (declination) 
due to drift-scan mode of observation each OGLE-II field is 
divided into 64 subframes before processing. Subframes are 
2048x128 pixels with a 14 pixel margin on each side. 

We compute the pixel posi tions of stars in the im- 
ages using the DoPHOT package l|Schechter. Mateo fc Salial 
1993). At the start of the processing for each exposure, the 
positions of stars in a single subframe are measured and 
cross-referenced with those in the template and the overall 
frame shift is obtained. Using this crude shift we can iden- 
tify the same region of the sky (corresponding to a given 
subframe of the template) throughout the entire sequence 
of frames. 

To treat properly spatial PSF variations, each 2048 x 128 
pixel subframe is divided into 4 smaller chunks with a size of 
512 x 128 pixels with a 14 pixel margin on each side. Then the 
positions of stars in all chunks are computed by DoPHOT. 
We use all stars with I < 18 mag (~ 400 of them, depending 
on the stellar density in each field) categorized by DoPHOT 
as isolated stars (marked as type=l) in the following analy- 
sis. We do not use the data points categorized by DoPHOT 
as a star blended with other stars (marked as type=3). 

The stars in each of the chunks are combined into the 
original 2048x128 subframes. We cross-reference the stars 
in the template and other frames with a search radius of 
0.5 pixels and derive the local transformation between these 
pixel coordinate systems for each subframe. We use a first 
order polynomial to fit the transformation. The resulting 
piece-wise transformation adequately converts pixel posi- 
tions to the reference frame of the template. Typical resid- 
uals are at the level of 0.08 pixels for bright stars (/ < 16) 
and 0.2 pixels for all stars (/ < 18). 

By using these transformation matrices, we cross- 
reference the stars in the template and other frames with 
a search radius of 1.0 pixels instead of 0.5 pixels used in 
ISumi. Ever fc Wozniakl i2003l) to increase the range of de- 
tectable high proper motion objects. We estimate that the 
probability of the mis-identification in this search radius is 
negligible (0.26 %). 

We have found that there are systematic differences in 
the mean positional shifts of stars from the template posi- 
tion {dx) and {dy) depending on time and pixel coordinate 
in x. We have measured the {dx) and {dy) of the stars in 81 
strips {X = ~ 80) centered at equal intervals in x coor- 
dinate between < x < 2048 with a width of ±25 pixels. 
Each strip contains typically ~ 2, 000 stars. In the upper 
panel of Fig. we show {dx) as a function of time for the 
strip X = 40 (x = 1024±25 pixels) in BULJ3C2. We can see 
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Table 1. Center positions of OGLE-II 49 fields. The number of frames Nf and stars JV S , and the mean of uncertainty in proper motion, 
< > (masyr - 1 ), are given as a function of / (mag) for each field, cr^ is averaged over I ± 0.5 mag. 
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the big jump at JD=2451041 (indicated by a vertical dashed 
line) where the exposure time of OGLE-II in the GB fields 
has been changed from 87 sec to 99 sec in the middle of 1998 
season, on August 15. In the upper panel of Fig.|2]we show 
typical mean positional shifts in x, {dx} of stars in strips 
in BULJ3C2 as a function of pixel coordinate x. The filled 
and open circles represent the (dx) of the frame taken at 
JD=2450887.822 (before the jump) and 2451336.769 (after 
the jump), respectively. There are also systematics in (dy) 
with the level of 0.04 pixels. We cannot see any such sys- 



tematics as a function of y pixel coordinate. Because of the 
good coincidence between the jump and the change in the 
drift scan rate that determines the effective exposure time, 
the bulk of the systematics may be caused by the change in 
drift scan rate. However, we do not know the detail reasons 
behind this at the present time. 

Even within the period before and after the jump, the 
shapes of Fig. differ from time to time and from field to 
field at the level of 0.04 pixels. By interpolating these curves 
of (da;) and (dy) as a function of x for each frame (time) of 
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Figure 1. Upper panel: The mean positional shift in x, (dx) of 
stars in a strip X = 40 (x = 1024 ± 25 pixels) in BUL.SC2 
as a function of time. We can see the big jump at JD=2451041 
(indicated by a vertical dashed line) where the exposure time of 
OGLE-II in the GB fields has changed from 87 sec to 99 sec in 
the middle of 1998 season, on August 15. Lower panel: The same 
figure as above after the systematic correction. 



■SC5 J l-6 ID=5935 1= 141093 V-I = '1.260 N=392 1 
-&i o „// J ) = (l-194(0.405), 0.939(0.346)) (mas/yr), a=-22.418 (mas) 
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Figure 3. Time variation of the position in acos<5 (upper) and 
<5 (lower) for star ID=5935 (V-I=1.260) in BUL.SC5. Filled 
circles represent actual positions, while open circles are posi- 
tions corrected for differential refraction with an offset of —0.16 
pixels. Solid lines indicate a model fit for the proper motion 
(jl a »,Hg) = (-1.19(0.41), 0.94(0.35)) (masyr" 1 ) with the 1 a 
errors in bracket. The differential refraction coefficient is a = 
-22.42 mas. 
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Figure 2. Upper panel: Typical mean positional shifts in x, (dx) 
of stars in strips in BUL_SC2 as a function of pixel coordinate x. 
The filled and open circles represent the (dx) of the frame taken 
at JD=2450887.822 (before the jump) and 2451336.769 (after the 
jump), respectively. Lower panel: The same figure as above after 
the systematic correction. Note that vertical scale in the lower 
panel is very different than in the upper panel. 



each field, we correct dx and dy for each star and frame. In 
the lower panel of Fig. and Fig.|5J we show the same plots 
after this systematic correction. This procedure is based on 
the assumption that average proper motions of a large num- 
ber of stars in separate groups of columns (i.e. different val- 
ues of X) should be the same. Note that the integral of the 
curves shown in Fig. [5] over all x-columns is unity, as this 
corresponds to the average position of all stars. 



An example of time dependence of the position for a 
star with a moderate detectable proper motion is shown in 
Fig.dlwith filled circles. To measure the proper motions in 
right ascension (fi a * = HaCosS) and in declination (^s), we 
fit the positions as a function of time t with the following 
formula: 



a = ao + Hu*t + a sin C tan z, 
S = So + fist + a cos Ctan z, 



(1) 
(2) 



where a is the coefficient of differential refraction, z is the 
zenith angle, and C is the angle between the line joining the 
star and Zenith and the line joining the star with the South 
Pole, and ao and 5o are constants. The parameter a is a 
function of the apparent star colour. We neglect the paral- 
lactic motion due to the Earth's orbit because its effect is 
strongly degenerate with the effect of differentia l refractio n 
for stars in the direction of the GB fever fc Wozniak 1200 ll) . 

In Fig. |21 we present with solid lines and open circles 
the best fit model for the proper motion (fj, a *,^s) and po- 
sitions allowing for the differential refraction respectively. 
As is written in this figure, the parameters for this object 
are: field SC5-1-6, which means that this object is in the 
OGLE-II field BUL.SC5, and chunk (X chmk ,Y chmk )=(l, 
6), OGLE ID=5935, / = 14.093, V -I = 1.260, the num- 
ber of data points = 392, proper motion (fj, a * , fJ.s) = 
(1.19(0.41), 0.94(0.35)) (masyr -1 ) with 1 a errors in brack- 
ets. The differential refraction coefficient is a = —22.42 mas. 

We computed ao, So, /J, a *, and a for all stars used 
to transform coordinate systems (approximately the num- 
ber of fields times the number of chunks times the typical 
number of stars per chunk, i.e. 49 x 256 x 400). In cases 
where the star is measured in the overlap region of more 
than one chunk of a given field, the data set with the largest 
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Figure 4. Colour Magnitude Diagram of stars with (T M < 2.5 
masyr -1 in BUL_SC2. Iq and (V — I)o are extinction corrected 
/-band magnitude and V — I colour. The stars in GB are defined 
within the ellipse centered at the center of RCGs plus 0.4 mag in 
I. 



number of points is selected. Stars with the fewer than 20 
data points are rejected. The catalogue of whole BUL_SC1 
was reco mputed in this analysis becaus e the catalogue pre- 
sented bv lSumi. Ever & Wozniakl J2003) contains only 70 % 
of this field. 

Our catalogue contains 5,080,236 stars, which is 79.8 % 
of all objects with I < 1 8 mag in the original OGLE tem- 
plate dUdalski et al.l2002l) . 7% of these stars do not have any 
1^-band photmetry. The missing V-band photometry of stars 
can be estimated from the differential refraction coefficient 
because there is good correlation between the differential re- 
fraction coefficient a and the apparent V — I colour for stars 
(see Fig. provided they belong to the same population 
as that of the majority. 

We have measured the mean proper motions of stars in 
the Galactic bulge (GB) defined in the ellipse in the Colour 
Magnitude Diagram (CMD) in Fig. 2] where the extinction 
a nd reddening are corrected by using the extinction map 
of ISumil (120031) . T he ellipse is l ocated at the center of the 
RCGs estimated in lSumi I fcOOot plus 0.4 mag in I and have 
semi-major and minor axises of 0.9 mag and 0.4 mag, respec- 
tively. This ellipse includes RCGs and Red Giants in the GB. 
Here we chose only objects whose proper motion accuracy 
is better than 2.5 masyr -1 . These mean proper motions of 
stars in the GB are assumed to be constant within each field. 
In Fig. we show these measured mean proper motions of 
stars in the GB as a function of x (upper panel) and after 
the systematic correction (lower panel) . We can see how well 
the systematic distortions are corrected. Note that our in- 
strumental reference frame is defined by all stars, and this 
is why the average proper motion of stars in the GB are not 
zero. 
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Figure 5. Upper panel: The mean proper motions in x (filled cir- 
cle) and in y (open circle), for stars in 50 pixels strips in BUL_SC2 
as a function of pixel coordinate x. Lower panel: The same figure 
as above after the systematic correction. 



4 HIGH PROPER MOTION OBJECTS 

The detectability of high proper motions is limited by the 
search radius used for cross - identification of stars on all 
images. The search radius, 1 pixyr - , corresponds to ~ 400 
masyr" 1 . Objects with n > 100 masyr -1 cannot be identi- 
fied over the full four year long observing interval of OGLE- 
II, as they move out of the search radius. In order to be 
able to follow fast moving stars over four observing seasons 
we made additional astrometry for objects for which pre- 
liminary estimates gave a proper motion fi > 100 masyr" 1 . 
Whenever a star moved more than 0.4 pixels from the pre- 
vious search center we moved that center to the median lo- 
cation of the last 3 data points. This procedure was adopted 
when it allowed us to locate the star in a larger number of 
CCD images. 

We show the positional movement of one of the high- 
est proper motion objects in Fig. [S] and images at 1997 and 
at 2000 in Fig. |7| This star has relatively fewer data points 
because this star is overexposed on the good seeing frames 
in the /-band. This star dose not have any V-band mea- 
surements because of the failure of cross-identification in 
the l/-band template image due to its high proper motion, 
/-band measurements of such objects may also be unreli- 
able because the OGLE photometric maps are based on 
the measurements b y the "fixed position mode" over years 
JUdalski et all2 002). The photometries of this star obtained 
by hand relative to the neighboring stars are: / = 11.70, 
V — I = 2.86. This colour is very red as expected from the 
relation between differential refraction coefficient a from the 
fit and colour in Fig|H| but its a and V — / do not match the 
relation exactly. 

The relation between a and V — / for BUL_SC42 in 
which the extinctions are relatively small, is slightly different 
from that for BUL_SC5 plotted in Fig. [H] The zero-point of 
a is 15 mas higher than that for BUL_SC5, but the expected 
colour of this star from this relation is still redder (V — I = 
4 — 5). In these figures the slopes differ between fields for 
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Figure 6. Same plots as Fig. |3] for ID=133638 (V-I=none) in 
OGLE-II field BUL_SC42, one of the highest proper motion stars. 
Filled circles represent the actual positions and open circles are 
the positions corrected for differential refraction, with an offset of 
+1.5 pixels. Solid lines indicate a model fit for the proper motion 
(/W>A<«) = (-261.19(2.51), -436.90(1.2)) (masyr- 1 ) with 1 a 
errors in bracket. The differential refraction coefficient is a = 
26.04 mas, which indicates the star is very red. This star dose 
not have y-band magnitude because of the high proper motion. 
The photometries of the star obtained by hand relatively to the 
neighboring stars are: / = 11.70, V — I = 2.86. 
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Figure 7. Images of the high proper motion star in Fig. i.e. 
ID=133638 in OGLE-II field BULJSC42, at 19 April 1997 (left) 
and at 30 September 2000 (right). 



the region V — I > 3. This is because of the difference of 
the population. In BUL_SC5, the majority in this colour 
region are RCGs and Red Giant Branch stars, but Red Super 
Giants in BUL_SC42. So this disagreement might be because 
this object is a nearby very red dwarf of M4-5 spectral type, 
not typical for this catalogue. 



5 CATALOGUE 

A sample for our catalogue of proper motions is shown 
in Table H The complete list of all 5,080,236 stars is 
available in electronic format via anonymous ftp from the 
server ftp://ftp.astrouw.edu.pl/ogle/ogle2/proper_motion/ 
and ftp://bulge.princeton.edu/ogle/ogle2/proper_motion/ 
The list contains star ID, the number p a * and (is and in 



galactic coordinate fii and fib with their errors, differential 
refraction coefficient a, standard deviation (Sdev) of data 
points in the fitting, equatorial coordinates «2uoo, £2000 in 
2000, apparent I-band magnitude and V — I colour, pixel 
coordinates on CCD x and y, and position of chunk X c and 
Y c to which this object belongs. ID, 2000 coor dinates, I and 

V — I for each object are identical to those in lUdalski et alJ 
(2002). If the object doesn't have V-band photometry the 

V — I colour is written as 9.999. 

Note that the number of data points TV differs from star 
to star even if they have similar brightness. This is because 
some are near the edge of the CCD image or CCD defects, 
and others are affected by blending. We use the positions 
which are categorized as a single isolated star by DoPHOT. 
Hence, blended stars can not be measured when the seeing 
is poor. 

In Fig. |S] we present a colour - magnitude diagram 
(CMD) for stars in OGLE field BUL_SC5, which is one of 
the most reddened OGLE-II fields. We also show the cor- 
relation between the differential refraction coefficient a and 
the apparent V — I colour for stars with I < 16 in this 
field. This correlation for other fields is similar, but it has 
a slight dependence on the population of the majority of 
stars in each field. In Fig. El we plot the uncertainties in fi a , 
(o>a„! upper panel), in fig (o> 5 , middle panel) and differ- 
ence between them (cr MQt — <r Mi5 , lower panel) in BUL_SC3 as 
a function of the apparent /-band magnitude. In the lower 
panel, where open circles represent mean values for each 1 
magnitude bin, we can see that <7 Mq „ is systematically larger 
than <r M5 at ~ 0.1 masyr -1 level. We see the same trend in 
all our fields. This trend is expected from the residual scat- 
tering due to the systematic correction and the differential 
refraction which is large in the direction of a. The mean 
uncertainties (<r M ) = (yo^ + af l5 ) with 2a clipping as a 
function of / are listed in Table Note: that the accuracy 
of proper motions in our catalogue is better than 1 masyr -1 
for 12 < I < 14. 

Fig. HOI shows the histogram of our proper motion mea- 
surements /x for the whole catalogue, with the lines with 
increasing thickness corresponding to all stars, and to those 
with the proper motion detected with confidence better than 
3a, 5a and lOer, respectively. The total number of stars with 
proper motions measured to 3a, 5a and 10cr accuracy is 
N 3 = 1, 469, 838, N 5 = 568, 665, N w = 61, 231, respectively. 
The inclined solid line: (log(N) — ~3log(/j,) + const.), has a 
slope corresponding to the expectation for a uniform distri- 
bution and kinematics of stars in space. The distribution of 
accurate (5a and 10cr) proper motions seems to be roughly 
consistent with the uniform distribution. 

However, the number of very high proper motion stars, 
(fj, > 200 masyr - ), appears to be smaller than expected 
from a uniform distribution. This may be due to saturation 
of images of nearby, and therefore apparently bright, stars 
in OGLE images. 

The distribution of less accurate (< 3a) proper motions 
has an apparent cut-off at n ~ 100 masyr -1 . Such stars are 
apparently faint, their positions have large errors, and they 
may be difficult to identify in a search radius of 1 pixel, 
which corresponds to fi ~ 100 masyr -1 . 

A more thorough analysis of various selection effects is 
beyond the scope of this paper. Readers are advised to use 
caution in a statistical analysis of our catalogue. 



Table 2: Sample of proper motion catalogue for BUL_SC2. 
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Figure 9. Uncertainties in ii a * (cr^ 



upper panel), fig (a^ 



panel) in OGLE-II field BUL_SC3 as a function of the /-band 
magnitude. In the lower panel, where open circles represent mean 
values for each 1 magnitude bin, we can see that cr^ at is system- 
atically larger than cr l _ ls by ~ 0.1 masyr level. 
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Figure 10. Histograms are shown for /i for all N = 5, 080, 236 
stars, and for measurements with confidence better than 3a, 5<r 
and IOct (from the thinnest to the thickest line, respectively) for 
all 49 OGLE-II fields. The total number of stars with proper 
motion with 3<r, 5<r and lOer accuracy are ./V3 = 1,469,838, 
N 5 = 568,665, Nio = 61,231. The inclined solid line: (log(N) = 
—3log(u) + const.) has the slope expected from stars that have 
uniform distribution and kinematics. 
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Figure 11. Mean proper motions of stars at GC (in) (filled circle) 
and (fif,) (open circle) within AI = 0.1 mag bin as a function of 
/-band magnitude. Solid and dashed lines represent the mean of 
these bins without weighting by errors. We can see the streaming 
motion of bright (I ~ 14.1) and faint (~ 14.7) RCGs. The proper 
motion of Red Giants (/ > 15) is the mid of them as expected 
from that they are on average on the GC. 



able in electronic format via anonymous ftp with the main 
catalogue (see §0. 

The proper motion measurements in our catalogue can 
be transformed to the inertial frame by formula; 



/iatOGLE — Ma* " Ma*,GC + Ma* ,GC, inert 
AWOGLE = flS — M<5,GC + M<5,GC, inert- 



(3) 
(4) 



Here (p a *, G c, inert, M<5,GC,incrt) = (-2.93,-5.17) masyr 1 is 
the expected proper motion of the GC relative to the inertial 
frame, assuming a flat rotation curve of v r ~ 220 kms -1 , 
the distance between the GC and the Sun of Rq = 8.0 kpc 
feisenhauer et al~l 12003') and Solar velocity of (v@i,VQb) = 
(5.25 kms" 1 , 7.17 kms -1 ) relative to the Local Standard of 
Rest (RSL) JPehnen fc Binnevl H998). This transformation 
gives us crude absolute proper motions, while this gives us 
very good relative zero points from field to field as we show 
later. 

The reader can apply this transformation formula to 
our catalogue to get value in the inertial frame. We didn't 
apply this transformation to our catalog because the trans- 
formation to the inertial frame can be improved with the 
QSOs to be discovered behind our fields in the future. 

To check our measurements we cross- identified stars i n 
our catalogue with the Tycho-2 catalogue jHog et alJl200Cl) . 
We selected from our catalogue objects with proper motions 
higher than 10 masyr" 1 and measured with a significance 
above 3<r, to avoid mis-identification. Most high proper mo- 
tion stars in Tycho-2 are saturated in OGLE images. We 
found 65 Tycho-2 stars in our catalogue and their proper 
motions: p a » (thin) and us (thick) are presented in Fig. 
1121 Here OGLE proper motions have been transformed into 
the inertial frame by equations and I^J. We can see a 
very good correlation between OGLE and Tycho-2 measure- 
ments. Dashed lines indicate hogle = HTycho-2, and solid 
lines represent best fit with fixed unit slope and a possible 
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Table 3. Proper motions of the GC for all 49 OGLE-II fields. The number of stars used in measurements N s tar, proper motions of the 
GC in our reference frame in equatorial coordinate li a ,GC an( i W,GC an d galactic coordinate /i;,GC an d Mb.GC with their errors. 



field 


Nstar 


Ma 


GC 


C7 ji t 




A^gc 


BULJ3C1 


23938 


-0. 


.57 





.07 


-0.55 


BUL.SC2 


25738 


-0 


.61 


0. 


.05 


-0.64 


BUL.SC3 


43597 


-0. 


.66 





.05 


-0.68 


BUL.SC4 


42847 


-0. 


.57 





.06 


-0.62 


BUL.SC5 


12704 


-0 


.72 


0. 


.08 


-0.94 


BUL.SC6 


10167 


-0 


.61 


0. 


.05 


-0.80 


BULJ3C7 


9431 


-0 


.55 


0. 


.06 


-0.68 


BUL.SC8 


7397 


-0. 


.98 





.05 


-0.97 


BUL.SC9 


7364 


-0. 


.85 





.05 


-1.02 


BUL.SC10 


8476 


-0 


.84 


0. 


.06 


-0.98 


BUL.SC11 


8067 


-0 


.80 





.06 


-1.03 


BUL.SC12 


9913 


-0. 


.95 





.05 


-1.33 


BUL.SC13 


9899 


-0. 


.84 





.06 


-1.13 


BUL.SC14 


19398 


-0 


.57 


0. 


.05 


-0.70 


BUL.SC15 


16454 


-0. 


.53 





.05 


-0.79 


BUL.SC16 


15801 


-0 


.71 





.06 


-1.00 


BUL.SC17 


16665 


-0 


.68 





.05 


-0.97 


BUL.SC18 


23064 


-0. 


.60 


0. 


.05 


-0.73 


BUL.SC19 


21441 


-0 


.60 


0. 


.05 


-0.77 


BUL.SC20 


31933 


-0. 


.42 





.06 


-0.61 


BUL.SC21 


30486 


-0 


.55 





.06 


-0.71 


BUL.SC22 


31278 


-0 


.57 


0. 


.05 


-0.60 


BUL.SC23 


27679 


-0 


.59 


0. 


.06 


-0.70 


BUL.SC24 


27355 


-0 


.55 





.06 


-0.59 


BUL.SC25 


25782 


-0. 


.52 





.05 


-0.66 


BUL.SC26 


21577 


-0 


.62 


0. 


.05 


-0.66 


BUL.SC27 


20250 


-0 


.60 


0. 


.04 


-0.72 


BUL.SC28 


12159 


-0. 


.73 





.05 


-0.78 


BUL.SC29 
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-0. 


.75 





.04 


-0.78 


BUL.SC30 


29140 


-0 


.54 


0. 


.05 


-0.64 


BUL.SC31 


28389 


-0 


.44 


0. 


.06 


-0.67 
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-0. 


.53 





.06 


-0.65 
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-0 


.61 





.06 


-0.69 
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-0 


.53 


0. 


.06 
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-0 


.56 


0. 


.05 
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.59 





.05 
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-0 


.68 


0. 


.06 
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.06 
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0. 
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-1.20 


BUL.SC48 


6965 


-1. 


.09 


0. 


.06 


-1.00 


BUL.SC49 


6299 


-0. 


.99 





.06 


-0.86 



offset. A good correspondence between Tycho-2 and OGLE- 
II proper motions gives us certain confidence in our mea- 
surements. Slightly larger offsets in /!# imply that the error 
in absolute proper motion is at a level of 1 masyr^ 1 , which 
can be improved by using QSOs in the near future. 

We also compared the measu rements in BUL_SC1 made 
bv ISumi. Ever fc WoznTakl 120031 (hereafter SCI') with the 
proper motions presented in this paper, for which our mea- 
surements are above the 10cr level of accuracy. The two sets 
of proper motions for 1368 cross-identified stars are shown 
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in Fig. H3l together with the best fit line. The offset between 
the two sets, and the differences between individual mea- 
surements are within estimated errors. Note: the scale is 
differen t than in Fig. 1121 Large reduce d chi square in [i a are 
because jSumi. Ever fc Woznia k 2003 didn't correct system- 
atic distortions (see though their systematic distortions 
have been reduced at a level of 2 masyr -1 by dividing im- 
ages into small chunks. 

We compared the measurements done by us in one of the 
overlap regions: between fields BUL_SC1 and BUL_SC45. 
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Figure 12. Comparison of proper motions fi a * (thin) and fig 
(thick) for 65 stars cross-identified in our catalogue and in Tycho- 
2 catalogue. The proper motions in our catalogue are transformed 
to the inertial frame by equations £3 and Dashed lines in- 
dicate fJ-oGLE = t 1 Tycho-2> an d solid lines represent best fit 
with the offsets: /i a *OGLB — fctTycho-2 = —0.23 masyr - 1 and 
WOGLE — WTycho-2 = 0.99 masyr -1 . fi a * are shifted for +30 
masyr - 1 for clarity. 

The proper motions of 115 cross-identified stars, and the 
best fit with a small offset between zero points, are presented 
in Fig. 1141 Here proper motions have been transformed into 
the inertial frame by equations J3J and Q. Good correla- 
tions between them with a rather small zero point offsets: 
Ma*sci - Ma*sc45 = 0.16 masyr -1 and fxssci — ossein = 
—0.07 masyr -1 gives us certain confidence in our measure- 
ments and in equations and J3) in terms of the relative 
the zero point. The scatter is also consistent with the esti- 
mated errors. 



7 POSSIBLE PROBLEMS 

There can be various problems associated with proper mo- 
tions of variable stars. Our fields are very crowded, hence 
many stars may be blended. In the case of blending we mea- 
sure an average position of a blend of several stars within 
a seeing disk. If a variable star is blended with other stars 
which have slightly different positions within a seeing disk, 
the average position of the blended image may change while 
one component of the blend varies. This change might mimic 
a proper motion. 

As an example we present the time variation of the po- 
sition of a ve ry long timescale microlensing event candidate 
JSmithll2003) in the top panel of Fig. [T51 The J-band light 
curve of this event is shown in the bo ttom panel of F ig. I15 | 
(ID=2859 in variable star catalogue of lWozniak et al . 2002). 
In the top panel we can see an apparent proper motion in the 
first year which is coincident with the apparent brightness 
fading, as shown in the bottom panel. After the second year 
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+100 -50 50 100 



fi scv (mas/yr) 

Figure 13. Same figure as Fig. ll2l for 1 368 stars cross-identified in 
BUL.SC1 bv lSumi. Ever fc Wozniakl2003l (hereafter SCI') and in 
SCI by this work, measured with better than 10<r level accuracy. 
Dashed lines indicate /^sci = Msd' which overlap with the solid 
line in fj,s, and solid lines represent the best fit with the offsets: 
Ha*SCl -PatSCV = °- n masyr" 1 and USSCX- USSCV = -0.03 
masyr -1 . ii a are shifted for +30 masyr -1 for clarity. Note: the 
scale is different from Fig. 1121 



the proper motion seems to be small, and the position data 
points are sparse, which coincides with the low brightness of 
the star and may reflect the presence of a blend - in poor 
seeing DoPHOT cannot resolve the two stars, hence there 
are very few data points in the upper panel. The blend was 
actually found in the higher resolution OGLE-III images in 
2002. 

A plausible interpretation is that the microlensed star 
was brighter than the nearby faint 'companion' star in 1997, 
but by 2000 it became the fainter of the two. The position 
of the 'companion' is consistent with the direction of the 
apparent proper motion. The faint stars seem to have low 
proper motion. High resolution observations are needed to 
fully understand this object. 

As another example we present in Fig. llul the time vari- 
ation of the position of a star ID=309705 in the OGLE-II 
field BUL_SC39. Filled circles represent measured positions 
with type=l (used in this work), which were at a fixed loca- 
tion for the first three years but moved significantly in the 
fourth season. The star ID=309705 identified these centroids 
at the edge of the search radius (y=-l) for the first 3 sea- 
sons and at the center (y=0) in the 4th season. On the other 
hand, the neighboring star ID=309653 which has a similar- 
brightness of I = 14.8 and position of 0.18 pixels East (pos- 
itive x) and -1.89 pixels South (negative y), identified the 
same centroids at the edge of the search radius (y=-l). 

To see the details of this object, in Fig. IKil we also plot 
the position measurements which are categorized as a star 
blended with other stars (type=3, which are not used in 
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Figure 14. Same figure as Fig. 1121 for 115 stars cross-identified 
in the overlap region between OGLE-II fields BUL_SC1 and 
BUL_SC45, and measured with better than 5<r level accuracy. 
Proper motions are transformed into the value in the inertial 
frame by equations J3J and J1J. Dashed lines indicate fJ-sci — 
fiSC45> an d solid lines represent the best fit with the offsets: 
l^a*SCl-^a*SC45 =0.16 masyr _1 and VSSCI- V6SC45 = -0.07 
masyr -1 . fj, a in SC45 are shifted for +30 masyr -1 for clarity. 
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this work) by DoPHOT for this star (crosses) and for a 
neighboring star ID=309653 (dots), which are shifted by 
+0.18 pixels in x and -1.89 pixels in y, i.e. these dots are 
as they are on the CCD, relative to ID=309705. We can see 
that these positions (crosses and dots) are identified around 
y=0 (ID=309705) and -2 (ID=309653), respectively, with 
filled circles in-between them. We also show the /-band light 
curves of ID=309705 (middle panel) and ID=309653 (bot- 
tom panel) in Fig. Iltil We can see ID=309705 is constant 
during four seasons but ID=309653 suddenly faded during 
the 4th season. This is likely to be a R CrB type variable. 

The simplest interpretation is as follows; small num- 
ber of data points in the seasons 1997-99 indicates that 
DoPHOT found that object which is composite of these 
two stars only on the bad seeing frames as a single star 
with type=l, in other cases they are separated but catego- 
rized as blended - type=3. In 2000 when the star ID=309653 
faded, the centroid of the composite moved and finally the 
star ID=309705 became a "single" star with type=l. So the 
number of data points is large in 2000. 

Proper motions of variable stars may have their errors 
increased not only because of variable contribution of blend- 
ing, but also because variable stars may change colours, and 
therefore the coefficient of differential refraction may also 
change. In rare cases of very long period variables this may 
be noticeable. Note: we do not treat this kind of objects in 
a special way, so a reader must be careful when using our 
catalogue in studies of variable stars. 

The effect of blending changes with seeing may con- 
tribute to the scatter of data points, but it is not likely to 



Figure 15. Top panel: Same plots as Fig. E] for the star 
ID=244353 in the OGLE-II field BUL.SC5. Filled circles repre- 
sent actual positions and open circles are correspond to positions 
corrected for differential refraction, with the offset of —1 pixels. 
Bottom panel: /- band light curve of same star (ID=2859 in cat- 
alogue of lWozniak et al.l200ll) . a possible very long microlensing 
event. 



have a seasonal effect. Hence, we think that seeing variations 
are not a major problem. 

The probability of blending is much larger for fainter 
stars, in particular those close to / = 18 mag. This may 
produce a bias in their proper motions, most likely reducing 
their formal proper motion, as the blended stars may have a 
different proper motion vector, so the average value is likely 
to be reduced. We also do not provide a special treatment 
for this kind of objects, so a reader must be careful when 
using our catalogue for faint stars. 



8 DISCUSSION AND CONCLUSION 

We have measured proper motions for 5,080,236 stars in all 
49 OGLE-II GB fields, covering a range of -11° < I < 11° 
and —6° < b < 3°. Our catalogue contains objects with 
proper motions up to fi = 500 masyr -1 and /-band magni- 
tudes in the range 11 < / < 18. The accuracy of proper 
motions in our catalogue is better than 1 masyr -1 for 
12 < / < 14. 
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One should keep in mind that all measurements of /i 
presented here are not absolute, but relative to the astro- 
metric reference frame which is roughly that of the Galactic 
Center (GC) with a small offset seen in the lower panel of 
Fig. 03 However, as demonstrated in iJS] by using the crude 
estimation for the proper motion of the GC in our refer- 
ence frame and formula equations and Q, we can ob- 
tain crude proper motions in an inertial frame. From the 
comparison with these inertial values and the Tycho2 cat- 
alogue, this transformation seems to work well with errors 
at a level of 1 masyr^ 1 . From comparison of fi measured in 
the overlap region of fields BULJ3C1 and BUL.SC45 (Fig. 
I14H . this transformation works very well in the relative off- 
set from field to field. These zero points for proper motions 
can be improved by using background quasars which may 
be detecte d in the near future using the OGLE-II variability 
catal ogue llWozniak et al.l2002l ; lEver l2002l : lDobrzvcki et alJ 
120031) . 

As demonstrated bv lSumi. Ever fc Wozniakl i2003Tl , the 
proper motions based on OGLE-II data can be used to 
clearly detect the presence of a strong streaming motion 
(rotation) of stars in the Galactic bar. While the reference 
frame established from all stars is not well defined with re- 
spect to the inertial frame, the relative motions of groups of 
stars within a given field are well determined. 

Though our primary goal is to constrain the Galactic 
bar model with the future analysis of our catalogue, we pro- 
vide proper motions for all stars with / < 18 mag in all 
49 OGLE-II GB fields because this catalogue might be use- 
ful for a variety of projects. An analysis of the catalogue is 
beyond the scope of the present study. 
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Figure 16. Top: Same plots as Fig. for star ID=309705 in 
OGLE-II field BUL_SC39. Filled circles represent actual positions 
with type=l (used in this work), crosses correspond to position 
measurements with type=3 (not used) and dots indicate positions 
measurements of neighboring star ID=309653 with type=3 which 
are shifted by +0.18 pixels in x and -1.89 pixels in y, i.e. these 
dots are as they are on CCD, relative to ID=309705. This is a 
probable case of a blend with one component being a variable 
star. 



